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I. INTRODUCTION 

There has been considerable interest, in recent years, in the transport properties of 
various materials in confined geometries (see, e.g., Refs. [l]-[5]). In many such applica- 
tions, boundary scattering of carriers can be expected to play a significant role. However, 
when efforts are made to calculate transport quantities such as the thermal conductivity 
for small samples, the effect of boundaries is often handled in a fairly simplified fashion. 
Boundary scattering of phonons, for instance, is typically incorporated by the addition 
of a constant term to the inverse scattering time for the phonon. The purpose of this 
paper is to provide a concise expression for the effect of boundaries on transport in small 
samples which improves upon this basic approach, without invoking the full machinery 
of the Boltzmann equation. Expressions are obtained for axial transport in samples of 
circular and rectangular cross section. These results will be applicable to samples which 
are small enough that the carrier mean-free path is on the order of the sample dimensions 
but not so small that the carrier spectrum is substantially modified from the bulk. In 
other words, the sample dimensions will be assumed to be much greater than the carrier 
wavelength. 

The method employed is based upon one proposed by Flik and Tien ||^ for the calcu- 
lation of the size effect in thin films. It is also related to a work put forward by Herring]^ 
in which approximate expressions of the thermoelectric power for some geometries are 
derived. The Flik and Tien method assumes that, for a carrier of a given frequency, a 
characteristic bulk mean-free path, /, can be defined which is determined by the other 
scatterers present in the sample. The goal, then, is to examine how this bulk value of 



2 



/ is modified by the presence of boundaries in tlie sample. Tlie calculation utilizes the 
concept of the exchange length l^^ ^, which is defined as the average distance normal 
to a plane that a carrier travels after having been scattered within that plane. To be 
more explicit, we consider a carrier that has undergone a scattering event within a plane 
that is perpendicular to the direction of net transport which will be referred to from here 
on as the positive z direction. We now allow the carrier to propagate to the point of 
its next scattering event, which, in the bulk, is a distance I away. This propagation is 
assumed to proceed with equal likelihood in all directions, lex is then defined as the aver- 
age z-component of all possible such propagation vectors, where the average is performed 
over the hemisphere in the positive z direction. The bulk value of this quantity, /qo is 
determined by the expression 0, 

loo = — / dcp dOlsinO cos9 = -, (1) 
271 Jo Jo 2 

where 6 is the angle between the vector / and the z axis and </> is the angle between the 
projection of / into the xy plane and the positive x axis. 

In the following, we consider the scenario in which the mean-free path is on the order 
of the sample dimensions. For this case, some carriers will strike the boundaries before 
traveling a full distance I and the exchange length will be correspondingly shorter. We will 
assume that scattering at the boundaries is diffuse, which will be valid when the carrier 
wavelength is smaller than the characteristic roughness features of the sample surface. 
We will also consider the sample to be free of grain boundaries, though this calculation 
could perhaps also be applied to samples whose grains have characteristic geometries 
which match those investigated here. A geometric analysis will be used to calculate the 
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average exchange length for axial transport in narrow samples (length is assumed infinite) 
of both circular and rectangular cross section. This quantity will then be related to an 
effective mean-free path, l^jj, which can be used to calculate axial transport coefficients. 
The calculations presented in the following sections are somewhat lengthy, but we assure 
the reader that the method pursued, particularly for the rectangular case, is considerably 
shorter than an approach one might naively take. 

II. CALCULATION OF THE EXCHANGE LENGTH 
A. Circulair cross section 

In order to calculate the exchange length for a cylindrical sample, we will initially 
assume that the excitation can originate with equal likelihood anywhere within a given 
circular cross section of the sample. The average value of the exchange length in the 
sample, l^x, will then be obtained by averaging l^.^ (which is itself an average over a 
hemisphere of solid angle) over the entire cross section. The geometry to be considered 
is shown in Fig. 1. We consider an excitation originating at some point a distance p 
from the center of the cross section of radius R and propagating in a random direction 
within the hemisphere of solid angle whose base is normal to the positive z direction. 
The quantity 9 is defined as the angle between the propagation vector, I and the z-axis, 
and (f) is the angle between the radius along which the origination point is located and 
the projection of I into the plane of the cross section. Note that I may or may not have 
length /, depending upon whether or not it is truncated by a boundary. We will also 
find it useful to define the quantity s{p, 0) as the length of the segment within the cross 
section plane which begins at the origination point and extends to the sample boundary 
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in the direction of (see Fig. 1). This will be the length of the projection of I into the 
plane of the cross section for the case in which the carrier hits the wall before traveUing a 
full mean-free path, The average exchange length will then be given by the expression, 

2 rR fir rTT/2 

/ea:=;^y^ c?PP dO sm 9 l^{p, (J), 9) , (2) 

where Iz is the z component of the propagation vector heading in the {9, 4>) direction 
and (f) has only been integrated over half its range for symmetry reasons. In the above, 
an integration over the polar angle within the cross section plane has already been 
performed since Iz does not depend on it. 

We now need to consider the separate cases in which I does and does not hit the wall. 
For a given p and 4>, I will hit within the range 9* < 9 < tt/2 where 9*, is given by 

r = sin- . 

The corresponding expressions for Iz for these two cases are given by 

/ cos ^ if < ^ < ^* 



h(p, 0, ^) = < 



(3) 

s{p,(j)) cot 9 if 9* <9 < 7r/2. 



We will also encounter situations, when I < 2R, where, for certain ranges of p and 
(f), s{p, 4>) > I. In these cases, I never hits the wall for any 9. How these ranges are 
determined will be examined more fully below. 

We are now in a position to evaluate the 9 integral. We will distinguish between the 
two cases described above, i.e. I sometimes hits the wall {Iz is described by Eq. (3)) and, 
I newer hits the wall {Iz — I cos 9 for all 9). These will be referred to hereafter as the SH 
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and NH regimes. The result of the 9 integration is 

f7r/2 



r desmeijp,(t),e) ^ < 

Jo 



1/2 never hits (NH), 

(4) 

s(p, 0) — s(p, 4>y/2l sometimes hits (SH). 



To proceed further, we must consider the relative magnitudes of I and R. There are 
two relevant possibilities, / > 2R and / < 2R. We will consider the simpler of these two 
cases first, I > 2R. In this regime, one can see that, with regards to the 9 integration, we 
will always "sometimes hit" the wall, for all possible values of p and (p. After performing 
the 9 integration, the expression for Ig^ becomes 

where, as can be determined from the geometry of the problem, 

s{p,(f)) = -p2sin2(/))V2 _pcos0. (6) 

It is possible to integrate the —8^/21 term in the order written above (0 first), but it 
can be shown that the s term must be evaluated by performing the p integration first to 
achieve an analytical solution. These operations yield the following result: 

Q p p2 

The case where I < 2R is somewhat more involved. Here we must distinguish between 
those regions of (p, (p) space where we are in the SH regime with regards to the 9 integra- 
tion and those where we are in the NH regime. The following results will be derived upon 
the assumption that I < R. The remaining case, R < I < 2R, can be shown, though it 
is not readily apparent, to yield results identical to the / < R case and we will therefore 



spare the reader their derivation. Assuming then, that / < R, one can see that in the 
region where p < R — I, I never hits the wall, regardless of (f) and 9. For p > R — I, on 
the other hand, we can be either in the SH or NH regime, depending upon the value of p 
and (f). It can be seen that, given a p greater than R — I, the angle 0* for which we cross 
over from the SH to the NH regime is determined by the condition that s(p, 0*) — I. The 
geometry for this situation is illustrated in Fig. 2. Using this condition, it can be shown 
that (f)* — TT — cos^^(^-^^^^-^). Summarizing these considerations for the I < R case, we 
arrive at the expression. 
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dpp [ dcl>-\. (8) 



The NH terms, where the integrand is 1/2, as well as the term involving — s^/2/, 
can be evaluated as presented above. However, the "s term" , as given by 

"s term" = ^ T dpp T c^0s(p,0), (9) 

TTit^ J R-l Jo 

can not be evaluated in the integration order shown, since the (f) integral is elliptic. So, as 
in the I > 2R case, we must reverse the order of integration and approach the p integral 
first. This procedure is less straightforward than the previous case, however, due to the 
fact that the upper limit of the 4> integral in Eq. (9) is p dependent. 

We must reconsider the problem from the perspective of integrating over p for constant 
(p. The relevant geometry for this case is presented in Figs. 3a and b. The goal is to 
determine the regions of phase space where we are in the SH regime, as this is the 
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only region which involves the s term that we seek to evaluate. If we first examine the 
different ranges of 0, we can see that for large values of this angle, we will never hit 
the boundary for any p or 6. The critical angle 02 above which this will be the case is 
determined by the condition that s{R, = I, (see Fig. 3a) and can be shown to be given 
hy (f)l — n — cos~^{^). Now, directing our attention to the different possible ranges of p, 
one can see that for < 02) we will be in the NH regime for small p and will cross over 
to the SH regime at a critical value p*. The determining condition for this crossover will 
again be that s{p*, 0) = /, (see Fig. 3b), which leads to the result 



These various considerations yield the proper expression for the s term with the integra- 
tion order reversed. 



where 02 and p* are given above. This expression can now be evaluated (see Appendix 
A) and, when combined with the results for the other terms obtained from Eq. (8), gives 
the total result for the I < 2R case: 



It is interesting to note that the —3^/21 term can not be integrated in a straightforward 
fashion with the p and integrations reversed, so both perspectives are necessary for a 
complete solution to the problem. 



p* = (i?2 _ ;2 _ I 




(10) 




I < 2R 



(11) 
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B. Rectangular Cross Section 

We now turn our attention to samples of rectangular cross section. The general 
expression for l^^ is quite analogous to the circular case and is given by 

2 /""/^ /■''/2 f^TT /•7r/2 

lex^ — r/ dy d(f) d9sm9 lz{x,y,(j),9), (12) 

nab Jq Jq Jo Jo 

where a and b are the lengths of the two sides and the integration is performed over 
only the bottom-left quarter of the cross section for symmetry reasons. The bottom-left 
corner is chosen as the origin of x and y and the a side lies along the x direction. We 
will define s{x,y,(l)) in a manner analogous to the previous case i.e., the distance from 
the origination point to the wall in the 4> direction, and point out that the results of the 
9 integration can again be divided into SH and NH regimes with the results shown in 
Eq. (4) applicable here as well. The expression for s will depend upon which boundary 
we are relating it to and, for the bottom wall for instance, is given by 

s{x,y,(j)) = y sec (j). (13) 

The integrations over 4>, x, and y require careful consideration to determine the regions 
of (x, y, (f)) space where we are in the SH or NH regimes. We will find that, as before, 
the makeup of these regions depends on the relative magnitudes of /, a, and h. We first 
examine the case when I > \/a? -\- W'. This is the simplest scenario because we are in the 
SH regime for all x, y and (f). The geometry for this case is illustrated in Figs. 4 and 
5. For a given origination point {x,y), we have to consider four ranges of 4> integration: 
01 — 02, 03 — ^ 04, 05 — ^ 06; 07 08, Corresponding to hitting the bottom, left, 
top and right walls, respectively, (see Fig. 4). For the bottom wall, is measured with 
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respect to the vertical through the origination point, with 0i in the negative direction 
and 02 in the positive direction. The counter-clockwise directional convention is chosen 
because I is directed into the plane of the paper. To avoid ambiguity in the evaluation of 
inverse tangents, which appear frequently in the calculation, a new zero of (p is defined for 
each side, with, for example, the zero with respect to the left side taken as the horizontal 
to the left of the origination point. 

We are now prepared to undertake the integration over the full range of cf). This task 
can be considerably simplified, however, by taking into account the symmetry of the 
problem. This process is illustrated schematically in Fig. 5. We begin by treating the 
contributions from the top and bottom, hereafter referred to as the "top/bottom" terms, 
separately from the sides. We will find that the latter can be related to the former by 
symmetry. Next, we note that the integration over the bottom-left quarter of our cross 
section of those 0's which will place us in the region where scattering will occur from the 
top boundary, i.e., 05 — > 06, is equivalent, by a 180 degree rotation, to an integration 
over the top-right quarter of the cross section of those 0's where scattering will occur 
from the bottom of the sample, i.e., 4>i —>■ 4>2- The part of our integral which accounts 
for scattering from the top and bottom boundaries then becomes 




The next simplification occurs as a result of an additional symmetry. It can be 
shown that when the integral is evaluated at its lower limit 0i, one obtains the same 
expression as when it is evaluated at the upper limit, 02, provided that the 0i expression is 
subjected to a change of variables which changes the region of integration in {x, y) space. 
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Specifically, the appropriate region is obtained by taking the original region and reflecting 
it about the line x — a/2. The benefit of this operation is that we reduce an expression 
which involves two terms, the integral evaluated at its upper and lower limits, to one 
involving only one term, simply by changing the region of (x, y) integration. In particular, 
we arrive at an expression which is the integral evaluated at its upper limit 02, where 
the region of (x, integration is the one given in Eq. (14) plus its reflection about the 
line X — a/2, (see Fig. 5). This leads to the following expression for the top/bottom 
terms: 



where ^2(2^, y) — tan~^(a;/y) and s{x, y, 0) is given by Eq. (13). One can go through the 
same logic for the scattering contribution from the sides and it can be seen, by making 
the substitution a; — > y and y — > x, that the expression obtained is identical to the one 
above with the exception that the upper limits of the x and y integrals are reversed. 
Summing these two contributions and performing the integrations we arrive at the total 
expression for the case in which I > y/a^ + 6^, 



The symmetry considerations illustrated here are not essential for performing the 
calculation above, but they become increasingly useful as one examines more complicated 
cases. Specifically, let us consider the case when I < a,b. As in the circular calculation, 
we will derive the results for a limited portion of this range and then demonstrate that 




(15) 




(16) 
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they are applicable for the entire range. We will initially require / to be less than a/2 and 
6/2. Taking, once again, the region of {x,y) integration to be the lower-left quarter of 
the cross section, one can see that / can only hit the bottom and left sides. The task is, 
as before, to discern which regions of {x, y, (p) space place us in the SH regime and which 
place us in the NH regime. The simplest way of visualizing this is to consider a circle of 
radius /, with its center located at various points {x, y) within the lower-left quarter, and 
to examine the manner in which it intersects the boundaries. We shall be particularly 
interested in the bottom wall, as the sides can be related to it by symmetry arguments. 
The relevant regions are shown in Fig. 6. In Region I, within a radius I from the bottom 
left corner, our circle intersects the bottom wall exactly once. For points within this 
region, when one considers the integration of over the bottom wall {i.e. from 0i to ^2), 
it can be seen that we are in the NH regime from 0i to some angle 0^ and in the SH 
regime from to 02- 0i is determined by the condition that s(a;, |/, 0^) = / and is given 



For points in region H, which have y < I but are farther than I from the corner, our circle 
intersects the bottom wall in two places, (f)\ and 02- 0i is as defined above and —02 = 'Pi- 
In this region, we find as we integrate across the bottom wall, that we are in the SH 
regime for 0* < < 02 and in the NH regime elsewhere. Finally, in Region HI, where 
y > I, our circle does not intersect the bottom wall at all and the NH expression applies 
to the entire integration. 

Now let us consider the top wall. It is clear that the NH expression applies as we 



by 




(17) 
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integrate from 05 to (pQ for origination points anywhere in the lower left quarter of the cross 
section. We can now apply the symmetry manipulations described for the I > y/o^~+W 
case to the present scenario. First, we rotate the "top" integration (05 to ^e) by 180 
degrees to convert it to an integration over the upper-right quarter from 0i to 02- We 
can now combine this contribution with the integral obtained for region III, as they 
both represent integrations from 0i to 02 where the NH expression holds for the entire 
range. The various integrals we need to evaluate, then, are shown in Fig. 7a, with their 
respective regions of {x, y) integration represented pictorially. These integrals represent 
the total contribution from integrations over the top and bottom walls to our expression 
for If^x-i i-e., the top/bottom terms. At this stage we can divide up the integrals into 
their constituent upper and lower limit parts, and by performing the type of variable 
substitutions and accompanying {x, y) integration region reflections discussed for the 
previous case (see Appendix B), reduce the problem to that shown schematically in 
Fig. 7b. The top/bottom contributions derived from Fig.7a are indicated by the diagrams 
on the left. Those on the right represent the contributions from the sides. To illustrate 
how these schematics are converted to mathematical expressions, we will write out the 
top/bottom terms (those indicated by the diagrams on the left in Fig. 7b) exphcitly: 

s{x,y,(l)f^ 



top/bottom — j dx j dy J d(l) is{x,y, (f)) — 



21 



=tan ^(x/y) 



+ 



+ 



dy 



.dx 



d(j) \^s{x,y,(f)) 



s{x,y,(t)f _ l_ 
21 2 



pa ro rl p\ I —x r ^ 

dx dy — dx dy / dcp- 

Jo Jo Jo Jo J 2 



(18) 



=tan ^{x/y) 



where, again, s{x,y,(l)) is given by Eq. (13). An analogous expression describes the side 
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wall contributions. The only changes are in the limits of the {x, y) integrations, and these 
are indicated by the diagrams on the right in Fig. 7b. Combining all contributions and 

performing the integrations it can be shown that the total expression for the case where 
/ < a, 6 is 

~ / ia + h)f /3 

Recall that the results above were derived for the case in which I < a/2,6/2. The 
case where b/2 < I < b appears, on the surface, to be quite different, since, in an analysis 
of the type described in Fig. 6, one would need to consider the possibility of hitting the 
top wall. However, an examination of the forms of the integration regions depicted in 
Fig. 7b reveals that they undergo no meaningful change until / becomes greater than b. 
This is one of the benefits of the manipulations which result in Fig. 7b; they reveal the 
fundamental equivalence of apparently disparate cases. The other principle benefit of 
the technique is that it obviates the need to extensively evaluate the various SH and NH 
regions for other, non-identical cases since the diagrams shown in Fig. 7b can be logically 
extended to other regimes. The case where 6 < / < a for instance, is shown in Fig. 8a. We 
note that the difference between the top/bottom diagrams and the side diagrams for this 
case results in an assymetrical dependence of lex on a and b. Performing the integrations 
indicated in the figure, we find that 

lex^F + G{b), (20) 



where. 
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COS 

TT 



\l) 3nql 



(21) 



and F is as defined in Eq. (19). 

The final case to consider is that for which a,b < I < \^a?~+lP. The appropriate 
diagrams are shown in Fig. 8b and they lead to the result, 



One final remark: the case where / > + 6^ which we analyzed first can now be seen 
to result from the final possible extension of the diagrams shown in Figs. 7 and 8. The 
shaded regions on the top line of Fig. 8b expand to encompass the entire rectangles and 
the terms on the second line go to zero. 



All of the above expressions, for both the circular and rectangular cases, were derived 
on the assumption of uniform origination, i.e., we assume that the excitation can originate 
anywhere within the cross section with equal likelihood. Recalling the physical interpre- 
tation of the exchange length, however, we note that the excitations we are considering 
are those that have just undergone a scattering event within the cross section plane. As 
pointed out in Ref. when the mean- free path of the excitation becomes much longer 
than the sample dimensions, the excitation becomes increasingly likely to scatter on a 
boundary and our assumption of uniform origination needs to be replaced with a bound- 
ary origination description. The expressions appropriate for boundary origination are as 
follows: 




(22) 



III. BOUNDARY ORIGINATION 
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Circular: 



2 rVa 
- / d<P 

7C Jo 
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T 



s{R,4>)- 



s{R, 



21 



(23) 



Rectangular: 



7r(a + b) [Jo 



+ I dy 
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dx 



tan ^{b/x) 







tan ^(a/y) 



is-s'/2l) + 



[s - s72/) - 

.tan-H^^) 



tan-l(^) 



d(t){s-s^/2l) 



-tan-i(2i) 



[s - s'/2l) 



7r{a + b) 



2 + b 



(a2 + 62) + 62 a6 



= J 



(24) 



7r(a + 6) TT tt/ 

The appropriate procedure is to match these solutions to those for uniform origi- 
nation at large Z's. This can be accomplished by a simple process which is as follows. If 
we wish the uniform solution, which we will call solution 1, to apply at some I = I' and 
we wish the boundary solution, which we will call solution 2, to apply for large /, the 
matched solution may be derived by the expression. 



matched solution = (sol 1) + exp(— m/'/Z)(sol 2 — sol 1). 



(25) 



where m is a matching parameter. We will take /' to be 2R for the circular case and 



Va2 + 62 for the rectangular case. The matched solutions then become 



Circular case: 



8R 
3^ 



— + exp(-m2i?//) 

2 ' M Stt 21 



I > 2R 



(26) 



16 



Rectangular case: 



■ex 



H + e^Y>{-md/l)[J - H] 



(27) 



where d = \/a?+lP, and H and J are defined by Eqs. (16) and (24), respectively. These 
forms result in small discontinuities at / = /' when they are combined with the solutions 
for / < /' derived in the previous sections, Eqs. (11) and (22). If the matching parameter 
m is chosen to be four, this discontinuity is less than 1.5%. It should be noted that, for 
the rectangular case when a ^ 6, we can expect to be largely in the boundary origination 
regime before I exceeds the length of the diagonal. The above equation does not allow 
for this possibility and is therefore most applicable to cases in which a is not too different 
from b. For samples where one dimension is much larger than the other, we refer the 
reader to the results of Ref. where an expression for thin films is derived. 



For the sake of clarity, the results for all the cases examined are summarized in Table 
I. Also included is the square subcase of the rectangular expression. 

I. Circular case: 



IV. DISCUSSION 



i.l<2R 




(28) 



a. I > 2R 




(29) 
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I. Rectangular case: 

i. I < a,b 

lex = F (30) 

a. b < I < a 

= F + G{b) (31) 

in. a,b < I < d 

L = F + G{a) + G{b) (32) 
iv. I > d 

lex = H + exp{-md/l)[I - H] (33) 

where F, G, H and / are defined in Eqs. (18), (20), (15), and (23), respectively. 

These expressions can now be used to aid in the calculation of transport quanti- 
ties in samples where boundary scattering is expected to play a role. We will examine 
the thermal transport in particular. In the kinetic theory approximation, the thermal 
conductivity is given by, 

K = -Gvl, (34) 
3 

where G is the contribution to the specific heat from the carrier in question and v is 
the carrier velocity. Now, in Ref. 0, it is asserted that, within this approximation, the 
transport Kz along the z axis in a sample of confined geometry can be obtained from the 
bulk thermal conductivity Koo and the exchange length along that axis by the expression 



It can easily be seen that this is equivalent to writing function of / as 

t^z{l) = f^ooihff) (36) 

with 

hff = 2/ea;- (37) 

In other words, the transport in the small sample can be obtained from the expression 
for the bulk transport by replacing / in the latter expression with l^jj as defined above. 

Because of the simple geometric nature of this argument, it is plausible that this sort 
of analysis can be applied to more sophisticated treatments of the thermal conductivity 
as well. For instance, in the case of phonon transport, the thermal conductivity is often 
written in the Debye approximation as an integral over phonon frequencies 0; 

where x is the reduced phonon frequency huo/ksT, ks is the Boltzmann constant, 6£, is 
the Debye temperature and r(T, x) is the frequency dependent scattering time. The total 
inverse scattering time, r(T, x)~^, is usually expressed as a summation of the inverse scat- 
tering times from scatterers of various types, e.g., point defects, phonon-phonon umklapp 
processes, etc. Within this context, the effect of boundaries is typically incorporated by 
adding a frequency independent term to this total of the form [|l^ 

T^^ = v/ad, (39) 

where d represents the sample dimension and a is a geometrical factor. 
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We propose that greater accuracy may be achieved from calculations involving Eq. (32) 
by omitting the boundary scattering term in the total inverse scattering time and utiliz- 
ing the above equations for the exchange length to modify the mean-free path instead. 
The proposed procedure is as follows: the factor r in Eq. (32) can readily be replaced by 
l{x,T)/v, at which point the integration can be seen to be over the frequency dependent 
mean-free path multiplied by another x dependent factor. For each such mean-free path, 
l{x,T), a value of lex can be derived by using the equation appropriate for the geometry 
of the particular sample under investigation. Each l{x, T) in the integral can then be 
replaced by an hffix, T) = 2lex{x, T) as indicated in Eq. (31). The proper, boundary lim- 
ited value of the thermal conductivity is then obtained by integrating over the hffix, T), 
with the other factors in Eq. (32) left unaltered. 

The dependence of lejj on / for various geometries is shown in Fig. 9. Also shown in 
Fig. 9 is the result when l^jj is determined by the expression 

1 I--' 



4i,= (^ + y) . (40) 



where a is taken to be 1 for the circular case and 1.12 for the square case [10|. This 
expression results from making the assumptions summarized in Eq. (33), i.e., the effect of 
boundary scattering is addressed by adding a constant term to the total inverse scattering 
time. One can see that the discrepancy is approximately 20% in the region where / is 
close to the sample dimensions. While this error is not terribly large for the square and 
circular cases, Eq. (34) can not be used at all for the rectangular geometry, where the 
sample can not be characterized by a single dimension. 

V. CONCLUSIONS 
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In summary, we have presented analytical expressions for the effects of boundary scat- 
tering in samples where the bulk carrier mean-free path is determined by other scatterers 
present. Results are derived for axial transport in long, narrow samples of circular and 
rectangular cross section, where scattering at the boundaries is diffuse. The results are 
incorporated into a definition of an effective mean-free path for axial transport which 
can be used to calculate coefficients such as the thermal conductivity. Though we have 
focused on thermal transport in the present work, the expressions derived here could be 
of use in the examination of a variety of transport phenomena in confined geometries. 
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APPENDIX A: Evaluation of Eq. (10) 

The evaluation of the integral expressed in Eq. (10) is not entirely straightforward 
and, for that reason, it will be discussed in greater detail here. The p integration is easily 
performed, resulting in the expression 



2 ^ 



"s term" = -— ( ' d(t> A{p, 0) 



7ri?2 Jo 

where 



■\/ B?—l'^ S\T? (f>—l COS (, 



(41) 



Aip, ^)-Jdpp[s--y - 3 Sin-/ - y ^''^ 
We have used Eq. (6) for s(p, 0) and written out the lower limit p* explicitly. The 
difficulty lies in the complicated nature of this lower limit. Simply substituting it into 
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A{p, 4>) leads to a result which is not easily dealt with. However, we can make use of 
other information to work this expression into a more tractable form. Specifically, we 
know that, from the definition of p*, when p = p*, s = I. Therefore, if A{p,(j)) can be 
expressed in terms of s, then, when the lower limit p* is substituted in for p, we can also 
substitute I for s. After some algebraic manipulations, A{p, 0) can be expressed as 

A{p, 0) = ( ^pI^SII + cos^ - pR^ cos - sys] . (43) 

sin y 3 J 

which, upon substitution of p* for p and I for s, becomes 

3 3 3 y sm y 

This can now be integrated with respect to 0, and, when combined with the upper limit 
term (which presents no special problems) gives the result 

APPENDIX B: Derivation of Fig. 7b 

In this appendix, we will elucidate the details of the process by which one goes from 
the form expressed in Fig. 7a to that shown in Fig. 7b. To aid in this discussion, a 
key of relevant diagrams, each representing a region of integration in the {x, y) plane, is 
provided in Fig. 10. We will also abbreviate the expression s — /2l diSW{s). 

We start by noting that the first term in square brackets on the first and second lines 
of Fig. 7a is the same. Therefore, in evaluating these particular two terms, their regions 
of (x, y) integration, /i and /2, can be combined to give region /y, (see Fig. 10). The 
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next step, as discussed in the main text, is to separate each integral into its upper and 
lower limit parts and, by performing a series of variable substitutions, reduce the problem 
to the smallest possible number of distinct terms. As an example of this dividing and 
substituting process, let us examine the term on the third line of Fig. 7a. This term can 
be written, with the help of Fig. 10, as 



where the multiplication notation indicates the double integration of the term in brackets 
over the {x,y) region defined as I3. Now consider the lower limit, 0i = — tan~^ (^^)" 
we make the change of variables x' — a — x, it can easily be seen that 



This variable change results, then, in a term that looks very much like the expression that 
one gets for the upper limit 02 = tan~^ with the exception that the region of {x,y) 
integration, Iq, is the mirror reflection of the original region, I^, about the line x — a/ 2. 
This procedure will now be apphed to the terms in Fig. 7a as follows: 





(46) 




(48) 



(47) 



(49) 



(50) 



(51) 



23 



In the above we have made use of the fact that 4>2 = —4>i, with 4>l given by Eq. (17). 
We note that one of the integrals had the change of variables applied to both limits (Eq. 
45) and one did not have it applied at all (Eq. 43) so as to generate the fewest number 
of similar terms. This number can be now be seen to be four: those involving W{s), 
evaluated at 02 and and those involving 1/2, evaluated at these same two limits. The 
final step in the derivation of Fig. 7b is the combining of terms, wherein the /'s of like 
expressions are patched together like pieces of a puzzle. This process is illustrated in 
Fig. 11 for the terms involving 1/2 evaluated at 02- It can be seen that four contributing 
diagrams, I2, I3, h, and Is, combine to yield Ig, which is the form that appears in Fig. 7b. 
The other top/bottom terms in Fig. 7b are derived in a similar manner. 
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Figure Captions 

Figure 1: Schematic showing the geometry relevant for the calculation of the exchange 
length for a wire of circular cross section. The case pictured is that for which I hits the 
boundary. Note however that the quantity s is always defined as the segment within the 
cross section plane that stretches from the origination point to the boundary, regardless 
of whether or not I actually hits the boundary. 

Figure 2: Schematic of the various regimes when / < R and R — I < p < R. The angle 
(j)* is defined by the condition that s(p, cf)*) — I. For angles greater than 0*, I never hits 
the boundary. Arrows indicate possible orientations of the vector I when 9 = n/2. 

Figure 3: a) Schematic indicating how the angle (p^, above which I never hits the boundary 
for any p, is determined. When (p — (p^, s{R,4>2) — I- b) Schematic indicating how p*, 
above which I sometimes hits the boundary, is determined for a given 0. When p — p* , 
s{p*,4>) = I. Arrows indicate possible orientations of the vector I when 9 = t:/2 for 
varying p and fixed 0. The relevant region for the p integration portion of the s term is 
indicated by the bracket. 

Figure 4: Diagram indicating the various values of which are used in the determination 
of the exchange length in a sample of rectangular cross section. 

Figure 5: Schematic representation of the symmetry processes used to simplify the cal- 
culation of the average exchange length. The rectangular diagrams represent double 
integrals in the (x, y) plane over the shaded regions. 

Figure 6: Diagram indicating the various regions in the (a;, y) plane which must be treated 
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in the calculation of the average exchange length. In region I, a circle of radius / intersects 
the bottom wall once. In region II, a circle of radius I intersects the bottom wall twice 
and in region III there is no intersection. The curve in the figure has radius /. 

Figure 7: a) Schematic representation of the various terms involved in the calculation of 
the contribution from the top and bottom walls to the exchange length for the case in 
which I < a/2,6/2. The rectangular diagrams represent integrations in the {x,y) plane 
over the shaded regions, b) Terms used to calculate the average exchange length lex 
after symmetry manipulations have been applied. The leftmost diagram in each term 
represents the contribution from the top and bottom walls and the rightmost diagram 
represents the contribution from the side walls. These terms must be multiplied by a 
pre-f actor to yield the final result for lex- The curves in the diagrams have radius I and 
are centered at the lower left corner. 

Figure 8: a) Terms used in the calculation of the average exchange length for the case 
in which b < I < a. The diagrams representing the {x, y) integrations can be seen to be 
logical extensions of the previous case. Again, the curves have radius /. b) Terms used in 
the calculation of the average exchange length for the case in which h,a <l < y/oF+W. 

Figure 9: a) Plot of the reduced effective mean- free path versus the reduced bulk mean- 
free path for the circular case. The solid line indicates the results when the exchange 
length method is used while the dashed hne indicates Casimir's results using Eq. (34) 
with a = 1. The inset shows a blow up of the low I region and has the same axes as the 
surrounding plot, b) Results for the rectangular case. The solid line shows the result for 
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a = b and the dot-dashed hne represents the case in which a = 2b. As in a), the dashed 
hne shows the result from Eq. (34) with a — 1.12. Inset is as in a). 

Figure 10: Key of diagrams used in Appendix B for the derivation of Fig. 7b. Curves in 
the diagrams have radius /. 

Figure 1 1 : lUustration of the way in which the various integration diagrams for the terms 
involving 1/2 evaluated at 02 combine to a single diagram. 
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Table Caption 

Table I: Summary of the equations for 1^^ for the various different geometries and the 
various ranges of the mean-free path /. R is the radius for the circular case, and a and 
b are the lengths of the sides for the rectangular case, with b taken to be the shorter of 
the two. d is equal to Vcfi + b"^- 



29 





R s^n crp 


f 

'•ex 


Circle 


/ < 2R 


8R B? j_ I „„„-! ( I \ J- r.r^o-l ( l'^ ^\ 

3^- IT + ^ COS (^iflj + l^cos [^-l) 

TT V 12 ' 24R2 ) 


I ^ zJrC 


8R B? I QU/1\ ( 4ii \ 

3^ - ^ + exp(-8ii://) y-^-^) 


Rectangle 


I < 


h (a, b) 


^ I < a 




(1 < 1 < (1 


r [a . (>) + Lj{n. 0) + Ct(o. « j 




if + exp(-4d//)(J-i/) 


Sim IP* "TP 


1 < n 




a <l < ^/2a 


F{a,a) + 2G{a,a) 


I > V^a 


f [2 + exp(-^)] [ln(l + V2) + (1 - V2)/3] 
-^[l + exp(-^)] 



Sirpq 127rpg 



(p,p2)3/2 



J 



7r(a+6) 



(t + In (^) + (1 + In (¥) 



+ 



27ri 



dr d ab 

7r(a+6) TT Tri 



30 



